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It is known that the thermal conductivity of a dilute gas can be derived by using 
kinetic theory. We present here a new derivation by starting with two known entropy 
production principles: the steepest entropy ascent (SEA) principle and the maximum 
entropy production (MEP) principle. A remarkable feature of the new derivation is 
that it does not require the specification of the existence of the temperature gradient. 
The known result is reproduced in a similar form. 
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INTRODUCTION 

Thermal conductivity, coefficient that appears in Fourier's law [H-Hl], 

3 = -«VT, (1) 
where j is the heat current density and T is the temperature. For a dilute gas, the thermal 

la, 

ncvX, (2) 



conductivity can be derived by using kinetic theory la in, with the result 
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where n is the number density, c is the molecular heat capacity, v is the molecular average 
velocity and A is the mean free path. This derivation is simple, but it requires one to 
specify the existence of the temperature gradient, VT. We know that, along the direction 
of VT, the average molecular energy increases at rate cVT. Also, for a given direction, a 
molecule needs to travel an average distance X z between two successive collisions. Thus one 
can simply consider along the VT direction a layer of dilute gas with thickness 2A 2 . Such 
a dilute gas can be viewed as two layers, each having height X z , and the two layers would 
exchange molecules. For each pair of molecules to be exchanged, the amount of energy that 
is exchanged is twice A z cVT. During unit time, across unit cross-sectional area, the number 
of molecules to be exchanged is nv z . Therefore the heat flux density is nv z \ z c^T . By using 
the relations \ z = iA and v z = -^=v, the thermal conductivity Eq. fl2]) is obtained. 

We shall present a new derivation. The new derivation does not require one to specify 
the temperature gradient. Instead, it uses the entropy production. Here, two aspects of 
the entropy production will be used: (1) for a given state, it evolves to a direction of 
the steepest entropy ascent; (2) among many candidates, the actual steady state has the 
maximum entropy production. In this paper, we call the first the steepest entropy ascent 
(SEA) principle and the second the maximum entropy production (MEP) principle. The 



name of "SEA" is originally introduced by Beretta [8l-lllj|. The MEP principle has been 
used by Paltridge fl2Tl4| to study the Earth's climate. There exist many studies of entropy 



production theory itself; see, for example [15|-|26|. This paper provides an application. 
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In order to derive the thermal conductivity of dilute gas, we shall consider a dilute gas that 
is fixed in between two thermal baths at different fixed temperatures. The dilute gas would 
carry a steady heat flux from whose expression the thermal conductivity can be extracted. 
We shall determine the form of the steady heat flux by using entropy production. In order 
to do that, first of all, we shall obtain the entropy of the dilute gas. 



ENTROPY OF DILUTE GAS 

In this section, we derive the entropy of a dilute gas. Let the dilute gas consist of N 
molecules, and let N be large. The entropy depends on the heat flux, which in turn is 
related to the microscopic configurations of each molecule, including in which direction it 
transfers energy, how much energy it transfers, and how often it transfers energy. 

Concerning the direction, we consider that each molecule transfers energy either upward 
or downward. Thus the probability that k molecules transfer energy upward obeys the bino- 
mial distribution B(N, |), which can be approximated by the normal distribution A/"(y, 

P(*) = ^exp(-^£!). (3) 

y^Y v 2 ' 

The height of the dilute gas is denoted by Az. We temporarily set Az equal to 2A Z so 
that the dilute gas can be viewed as two layers, each of which has height X z and consists of 
N/2 molecules. For the bottom layer, if k\ molecules transfer energy upward, the number 
of molecules that transfer energy downward would be N/2 — kx, while, for the top layer, k% 
molecules transfer energy upward and N/2 — k 2 molecules transfer energy downward. Here, 
k\ + &2 actually is k. Thus the energy transferred from the bottom layer to the top layer is 
proportional to 

h-(^-k 2 ) = k-^. (4) 

Concerning how often a molecule transfers energy, on average a molecule transfers energy 
once every time interval X/v, since on average a molecule collides once during that time 
interval. Concerning how much energy a molecule transfers, we denote the average of that 
by e. We know that, when two molecules collide, the energy transferred is a little bit less than 
the energy difference between them. The average of that energy difference is comparable 
with the molecular energy distribution width, which is further comparable with the average 
molecular energy cT. So e is comparable with cT, but a little bit less. 

Thus the heat flux is 

(5) 
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By using Eq. (j3J), we see that 



p(J)= 7^ exp pw J J' (6) 

This equation reflects the fact that, even for a dilute gas that is isolated, a fluctuating heat 
flux J may still occur with certain probability P(J)- 
We know that 

Q(J)ocP(J), (7) 
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where Q(J) is the number of the microscopic states that are associated with J. Using the 
Boltzmann formula 

S(J) = k B \nQ(J), (8) 
we write for the entropy of the dilute gas 

S(J) = S — J 2 = S — -— t^J 2 > (9) 

Ne z v 2 2 A ne 2 v 2 

where So is introduced to denote the entropy of the dilute gas when it is in the equilibrium 
state, N is equal to 2nAX z , and A is the cross-sectional area. Here, Az has been temporarily 
set equal to 2X Z . But notice that Az can be extended to any value, because entropy is an 
additive physical quantity. Thus we write the general form of the entropy of a dilute gas 
with any height Az as 

S(J) = S Q - \c 2 J\ (10) 

with 

C = *£%L. (11) 

A ne z v z 

We may also obtain the entropy per unit volume, s = t^-, with respect to heat flux 
density j = J /A, 

s(j) =so- i^j 2 . (12) 
2 ne z v z 

A similar formula is also proposed by Lebon et al. 27], |28( , in the form 

ps(T, q) = ps (T) - l -r{\T 2 )- l q ■ q, (13) 

where p is the mass density, s is the entropy per unit mass, q is the heat flux density 
vector, and r is the relaxation time; A has a different meaning here: it means the thermal 
conductivity. In the end, we will obtain the thermal conductivity ( 122|) . by using which one 
can see that Eqs. ( |T2l and ( ITBl are the same. 



ENTROPY PRODUCTION AND THERMAL CONDUCTIVITY OF DILUTE GAS 

When the dilute gas is fixed between the two baths at different fixed temperatures, a 
steady heat flux will flow through it. In order to determine the steady heat flux, it seems 
that, for every possible heat flux, we need to examine all possibilities of how it evolves. 
Fortunately, it will turn out that, if we use the SEA principle in conjunction with the MEP 
principle, for a given heat flux we only need to examine two specific possibilities. 

One possibility is that the given heat flux J is steady. For this possibility the entropy 
production would be 

o- = d - ±r)J, (14) 

J-l J-2 

where T\ and T2 are the fixed temperatures of the baths and T\ < T<i. Here, the heat flux 
flows from the bath at the higher temperature T2 to the bath at the lower temperature T\. 
The other possibility is that the given heat flux J starts to relax at the maximum rate, 

J R (t) = Jexp(-~), (15) 
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where r is the minimum possible relaxation time, which is approximately equal to X/v. To 
see why the relaxation is exponential and r ~ X/v, let us study an example. In the example, 
let the height of the dilute gas be 3A 2 , so that it can be viewed as three layers, each having 
height X z . The two end layers are in contact with the thermal baths and their temperatures 
are therefore different. Through all three layers a heat flux flows steadily. Now, suppose 
that the thermal baths are suddenly detached. The heat flux of the two end layers will 
cease immediately while the heat flux of the middle layer will relax gradually. Thus the 
two end layers act as new thermal baths but their temperatures are brought closer and 
closer by the heat flux flowing through the middle layer. The closer the two temperatures 
are, the smaller the heat flux becomes. The smaller the heat flux becomes, the slower the 
rate that the two temperatures get close will be. Thus both the temperature difference and 
the heat flux decrease exponentially. Concerning r, one needs to consider the microscopic 
mechanism of the heat flux. One may view the middle layer as a collection of individual 
molecules bouncing between the two end layers. On average, each bounce takes time X/v 
and transfers energy comparable with the molecular energy difference between the two end 
layers. Thus just a few bounces can bring the two temperatures close. So the relaxation 
time t is comparable with X/v but a little bit larger. A dilute gas having height larger than 
3A 2 can be viewed as many layers, each having height X z . Then the discussion would be 
the same, but each layer except for the two end layers plays two roles: it carries a relaxing 
heat flux and it acts as a thermal bath. Thus all middle layers are correlated, and such 
correlation makes the relaxation time larger. But we are only concerned with the instant at 
which the heat flux starts to relax. At that instant, the correlation has not yet taken effect, 
and thus, in the limit t — > 0, r is the same for a dilute gas having any height. 
Using Eq. (JTU|) . we see that 

S(t) = S Q - l -C 2 J R {tf = S - l -C 2 J 2 exp (-2^) , (16) 



and 



dS(t) C 2 2 / t\ 
a(t) = — = -J exp ( -2- ) , (17) 



which, in the limit £ — >- 0, leads to 



a = — J 2 . (18) 

T 

The reason why we use the limit t — ^ is that we only need to know whether or not the given 
heat flux would start to relax. The relaxation process itself is irrelevant to our problem. 
(Also, in the limit t — > 0, Jr(£) should be uniform in space and time.) 

There is one issue to be addressed: the entropy production within each bath can be 
neglected if the baths are chosen to be perfect. By 'perfect' we mean that their thermal 
conductivity k is very large and their C 2 is very small. When their k is very large, the part 
of their entropy production arising in the form of Eq. (jl4p can be neglected, because the 
temperature distribution within each bath would be uniform. When their C 2 is very small, 
the part of their entropy production arising in the form of Eq. ( I18p can be neglected too. 
So we do not need to consider the entropy production of the baths. 

The entropy productions of the two possibilities are plotted in Fig. [TJ They are equal 
when J = J a- We shall use first the SEA principle and then the MEP principle to determine 
that the steady heat flux must be J a- According to the SEA principle, for a heat flux J > J a, 
the heat flux cannot be steady because the entropy production would not be the maximum. 
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FIG. 1: Two entropy productions with respect to small heat flux. The solid line is the entropy 
production when the heat flux is steady. The dashed line is the entropy production when the heat 
flux starts to relax at the maximum rate. J a is the determined steady heat flux. For a heat flux 
J > Ja, it cannot be steady because in this range the solid line is not the highest, at least is 
lower than the dotted line. The steepest entropy ascent (SEA) principle indicates that a state can 
occur only if it is associated with the maximum entropy production. Thus the steady heat flux 
must be in the range < J < J a- Subsequently, according to the maximum entropy production 
(MEP) principle, we know that the steady heat flux must be Ja because among all candidates it 
is associated with the maximum entropy production. 



(But it would not start to relax at the maximum rate either. Actually, the heat flux will 
start to relax at a smaller rate until J = J a, and during the entire relaxation process the 
dilute gas will continue to exchange heat with the baths.) Thus the steady heat flux must 
be in the range < J < Ja- Then according to the MEP principle, we see that the steady 
heat flux must be Ja because, among all candidates, Ja is the one that is associated with 
the maximum entropy production. 

The steady heat flux in the figure is J a, and therefore the steady heat flux is 

J = (19) 

Plugging Eq. ffTTT) in, we write for the heat current density 

. _ J _ rne 2 v 2 11 1 

3 = A = ~ZhT A^ ( T\ ~ T 2 h ( 0) 
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which becomes, in the limits Az — > and T 2 — > Tx, 

rne 2 v 2 dT 

3= Zk^f 2 ~ltz~' ( ' 

From Eq. (12 ip and Eq. ( I14p . we see that the heat flux indeed flows from the bath at the 
high temperature to the bath at the low temperature. This is the same as the minus sign in 
the equation (CQ) indicated. Finally, the thermal conductivity of the dilute gas is extracted 
as 



2 2 

me v 



3k B T 2 



(22) 



If one uses the approximations 



r ~ — , e ~ cT, k B ~ c, (23) 



v 



one will find that the new result (122]) is the same as the kinetic theory result ([2]). In fact, 
e is a little bit less than cT, as discussed before Eq. ([5]), while r is a little bit larger than 
X/v, as discussed after Eq. (TT5]) . Also, c is a little bit larger than k B ', for example, for an 
ideal gas c = \k B . So the two results may be different, but they should remain comparable 
with each other. 



CONCLUSION 



We have derived the thermal conductivity of a dilute gas by using a new approach of 
entropy production. The derivation is based on the fact that there exists an entropy pro- 
duction competition between at least two specific possibilities. One possibility is that the 
heat flux is steady. The other possibility is that the heat flux relaxes exponentially at the 
maximum rate. Then, by applying the steepest entropy ascent (SEA) principle, we have 
obtained the upper limit of the steady heat flux. Subsequently, by applying the maximum 
entropy production (MEP) principle, we have identified the actual steady heat flux. The 
resulting expression for the thermal conductivity is comparable with the result of kinetic 
theory. As an intermediate step, we have found that the non-equilibrium entropy of a dilute 
gas is lower than the equilibrium value by an amount proportional to the square of the heat 
flux. 
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